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Abstract: We consider a class of point processes (pp), wliich we 
call sub-Poisson; these are pp that can be directionally-convexly 
(dcx) dominated by some Poisson pp. The dcx order has already 
been shown in [4] useful in comparing various point process 
characteristics, including Ripley's and correlation functions as 
well as shot-noise fields generated by pp, indicating in particular 
that smaller in the dcx order processes exhibit more regularity 
(less clustering, less voids) in the repartition of their points. 
Using these results, in this paper we study the impact of the dcx 
ordering of pp on the properties of two continuum percolation 
models, which have been proposed in the literature to address 
macroscopic connectivity properties of large wireless networks. 
As the first main result of this paper, we extend the classical 
result on the existence of phase transition in the percolation of 
the Gilbert's graph (called also the Boolean model), generated 
by a homogeneous Poisson pp, to the class of homogeneous sub- 
Poisson pp. We also extend a recent result of the same nature 
for the SINR graph, to sub-Poisson pp. Finally, as examples we 
show that the so-called perturbed lattices are sub-Poisson. More 
generally, perturbed lattices provide some spectrum of models 
that ranges from periodic grids, usually considered in cellular 
network context, to Poisson ad-hoc networks, and to various more 
clustered pp including some doubly stochastic Poisson ones. 

Index Terms — percolation, dcx order, Gilbert's graph. Boolean 
model, SINR graph, wireless network, Poisson point process, per- 
turbed lattice, determinantal point process, connectivity, capacity 

I. Introduction 

A network, in the simplest terms, is a collection of points 
in some space (e.g. on the Euclidean plane), called nodes or 
vertexes, and a collection of node pairs, called edges. The 
presence of an edge between two nodes indicates that they 
can directly communicate with each other. The mathematical 
name for this network model is graph. 

A class of networks that has recently attracted particular 
interest in the wireless conununication context, is called ad- 
hoc networks. It is distinguished by the fact that the network 
nodes are not subject to any regular (say periodic) geometric 
emplacement in the space but can be rather seen as a snapshot 
of some random point pattern, called also point process (pp) 
in the mathematical formaUsm typically used in this context. 

Connectivity, i.e. possibility of indirect, multi-hop, com- 
munication between distant nodes, is probably the first issue 
that has to be addressed when considering ad-hoc networks. 
An ubiquitous assumption when studying this problem is that 
the node randomness is modeled by a spatial Poisson pp. 
This latter situation can be characterized by independence 
and Poisson distribution of the number of nodes observed 
in disjoint subsets of the space. Poisson assumption in the 
above context is often too simplistic, however analysis or even 
modeUng of networks without this assumption is in most cases 
very difficult. 

In this paper we introduce some class of point processes that 
might be roughly described as exhibiting less variable point 
patterns than Poisson pp. We call point processes (pp's) of this 



class sub-Poisson point processes. Our objective is two-fold. 
On one hand we want to argue that this class of pp's allows to 
extend some classical results regarding network connectivity 
and on the other hand we want to bring attention to sub- 
Poisson pp's, as they might be useful for modeling of ad-hoc 
networks, whose nodes are more regularly distributed than a 
Poisson pp. 

A. Sub-Poisson Random Variables 

A random variable is called sub-Poisson, if its variance is 
not larger than its mean (with the equaUty holding true for 
Poisson variable). Intuitively, if strictly sub-Poisson variables 
were to describe the number of nodes in different subsets of the 
space then the resulting point patterns would exhibit less clus- 
tering (or bunching) than the Poisson point pattern having the 
same number of points per unit of space volume. This in turn, 
still intuitively, should have positive impact on the connectivity 
and perhaps capacity and other network performance metrics, 
the reason being that the perfectly regular, periodic patterns 
are commonly considered (and sometimes can be proved) as 
being optimal. The aim of the present article is to provide 
rigorous results on the comparison of the connectivity and 
capacity properties of certain sub-Poisson networks to these 
of the respective Poisson networks. 

B. dcx Sub-Poisson Point Processes 

The statistical variability of random variables (say with the 
same mean) can be compared only to some limited extent by 
looking at their variances, but more fully by convex ordering. 
Under this order one can compare the expected values of 
all convex functions of these variables. In multi-dimensions 
there is no one single notion of convexity. Besides different 
statistical variability of marginal distributions, two random 
vectors (think of number of nodes in different subsets of the 
space) can exhibit different dependence properties on their 
coordinates. The most evident example here is comparison of 
the vector composed of several copies of one random variable 
to a vector composed of independent copies sampled from 
the same distribution. Among several notions of "convex-like" 
ordering of random vectors the so called directionally convex 
{dcx) order, allowing one to compare expectations of all dcx 
functions (see Section II-A below) of these vectors, takes into 
account both the dependence structure of random vectors and 
the variability of their marginals. It can be naturally extended 
to random fields by comparison of all finite dimensional 
distributions, as well as to random pp' as and even locally finite 
random measures by viewing them as non-negative fields of 
measure-values on aU bounded Borel subsets of the space; 
cf. [4]. 

Using this latter formalism, we say that a pp is dcx-sub- 
Poisson (or simply sub-Poisson when there is no ambiguity) 



if it is dcx-smaller than a Poisson pp having the same mean 
measure, i.e., the mean number of nodes in any given set. We 
will also say sub-Poisson network in the case when the nodes 
of the graph modeling the network are distributed according 
to some sub-Poisson pp. 

We shall also see that there are classes of pp whose so- 
called joint intensities (densities of the higher order moment 
measures, when they exist; see Section II-B) are smaller than 
those of Poisson pp. Using this latter (and weaker) property 
we define the class of weakly sub-Poisson point processes. 

Another weakening of the comparability assumption, to 
the expectations of increasing (or decreasing) dcx functions 
allows to define increasing dcx sub-Poisson pp's, which can 
be seen as being less variable than some Poisson pp having 
possibly smaller (or larger) mean measure. We shall abbreviate 
increasing dcx by idcx and similarly decreasing dcx by ddcx. 

Our choice of the dcx order to define sub-Poisson pp's has 
its roots in [4], where one shows various results as well as 
examples indicating that the dcx order implies ordering of 
several clustering characteristics known in spatial statistics 
such as Ripley's K-function or second moment densities. 
Namely, a pp that is larger in the dcx order exhibits more 
clustering (while having the same mean number of points in 
any given set). 

C. Connectivity of Sub-Poisson Networks 

FuU connectivity (multi-hop communication between any 
two nodes) of an ad-hoc network with many nodes is typ- 
ically hard to maintain, and so, a more modest question of 
existence of a large enough, connected subset of nodes (called 
component) is studied. A possible approach to this problem, 
proposed in [10], and based on the mathematical theory of 
percolation, consists in studying existence of an infinite (called 
giant) component of the infinite graph modeling a network. 
Existence of such a component is interpreted as an indication 
that the connectivity of the modeled ad-hoc network scales 
well with its number of nodes. 

1} Percolation of Gilbert's Network: Percolation models 
have been extensively studied both in mathematical and com- 
munication literature. The model proposed in [10], called now 
the Gilbert's model \ is now considered as the classical 
continuum model in percolation theory. It assumes that each 
node has a given fixed range of communication p and direct 
connection between any two nodes is feasible if they are 
within the distance p from each other, regardless of the 
positions of other nodes in the network. The known answer 
to the percolation problem in this model is given under 
the assumption that the nodes are distributed according to 
a homogeneous Poisson pp having a density of A nodes on 
average per unit of volume (or surface). The result says that 
there exists a non-degenerate critical communication range 
< Pc = PcW < 00, such that for p < Pc there is no 
giant connected component of the network (the model does 
not percolate), while for p > p^ there is exactly one such 
component (the model percolates), both statements holding 
true almost surely; i.e., for almost all Poisson realizations of 
the network. An equivalent statement of the above result says 

'also Boolean model or random geometric graph 



that for a given communication range p there exists a non- 
degenerate critical density of nodes < Ac = Ac(p) < oo, 
below which the model does not percolate and above which 
it does so almost surely. 

As one of the main results of this paper we will prove 
an extension of the above result, which says that the crit- 
ical communication range of any homogeneous sub-Poisson 
Gilbert's network model is not degenerate. Moreover, this 
critical communication range is bounded away from zero and 
infinity by the constants which depend only on the mean 
density of nodes and not on the finer structure of the sub- 
Poisson pp of network nodes. Partial results, regarding only 
non-degeneracy at zero or at infinity can be proved for idcx— 
and ddcx-svto Poisson networks respectively. 

2) Percolation of SINK Networks: A more adequate per- 
colation model of a wireless communication network, called 
the SINK graph, was studied more recently in [6, 7]. It allows 
one to take into account the interference intrinsically related to 
wireless communications. The interference power is modeled 
by the shot-noise field generated by the pp of transmitting 
nodes. Each pair of nodes in the new model is joined by an 
edge when the signal power to this shot-noise plus some other 
(external) noise power ratio is large enough. The resulting 
random graph does not have the independence structure of 
the Gilbert's model and increasing the communication range 
(equivalent in this model to increasing the signal power and 
hence the value of the shot-noise) is not necessarily bene- 
ficial for connectivity. Similar observation holds as regards 
the increase in the node density. In fact, the above SINR 
network model has two essential parameters: the density A 
of Poisson pp of nodes and the shot-noise reduction factor 
7. The percolation domain is characterized in the Cartesian 
product of these two parameters. The key result of [7] says 
that whenever A is larger than the critical value corresponding 
to the percolation of the model with interference perfectly 
canceled out (for 7 = 0, which simplifies the model to the 
Gilbert's one) then there exists a critical value 7c = 7c(A) > 
such that the model percolates for 7 < 7c and does not 
percolate for 7 > 7^. As the second main result of this 
paper, we extend the above result, also to sub-Poisson SINR 
networks. 

D. Comparison of Shot-Noise Fields 

Our proofs rely on the comparison of extremal and additive 
shot-noise fields generated by dcx-ordered pp's studied in [4]. 
While the connection to the SINR model is evident (additive 
shot-noise is an element of this model) the connection to the 
Gilbert's model is perhaps less evident and relies on the fact 
that this latter model can be represented as a upper level-set 
of some extremal shot-noise field. 

More precisely, from [4, Propsition 4.1] one can conclude 
that the probabiUty of n given locations in the space not being 
within the communication range of any of the nodes in the 
Gilbert's model is higher for the network whose nodes are 
modeled by a pp larger in dcx order. Using this property, 
suitable discretization of the model and the Peierls argument 
(cf. [11, pp. 17-18] or [1, Proposition 14.1.4]) one can prove 
finiteness of the critical transmission range in the sub-Poisson 
Gilbert's network. The strict positivity of this range can be 



proved by comparing the expected number of paths from the 
origin to the boundary of an increasing box, again in some 
suitable discretization of the model. This latter comparison 
can be done relying only on weak sub-Poisson assumption. 

The remaining part of this paper is organized as follows: 
In Section II we provide necessary notions and notation. Sub- 
Poisson Gilbert's model is studied in Section 111 and sub- 
Poisson SINR model in Section IV. In Section V we make 
some remarks on the impact of dcx ordering on various model 
characteristics usually called capacities, including the so-called 
capacity functional, being one of the fundamental charac- 
teristics studied in stochastic geometry, as well as on some 
other capacity quantifiers in the information-theoretic sense. 
In Section VI we show some examples of sub-Poisson and 
weakly sub-Poisson pp, in particular the so-called perturbed 
lattices. Conclusions as well as open questions are presented 
in Section Vll. 

11. Preliminaries 

A. dcx Order 

We say that a function / : R'' — >^ R is directionally 
convex (dcx) if for every x,y,p,q such that p < x,y < q 
(i.e., p < X < q and p < y < q with inequalities understood 
component-wise) and x + y = p + q one has f{x) + f{y) < 
f{p) + f{q). Also, we shall abbreviate increasing and dcx 
functions by idcx and decreasing and dcx by ddcx. 

For two real-valued random vectors X and Y of the same 
dimension, X is said to be dcx smaller than Y, (denoted by 
X <dcx Y) if E(/(X)) < E(/(y)) for every dcx function 
/ for which both expectations are finite. In full analogy 
one defines idcx and ddcx orders of random vectors con- 
sidering idcx and ddcx functions, respectively. These orders 
clearly depend only on the distributions of the vectors. Two 
real valued stochastic processes (or random fields) are said 
dcx, idcx or ddcx ordered if any finite-dimensional distribu- 
tions of these processes are ordered. This definition extends 
also to locally finite random measures (in particular pp's), 
by ordering random values of these measures (in particular 
numbers of points) on any finite collection of subsets of the 
state space. More precisely, we say for two pp's $1, $2 that 
^"1 <dcx{idcx,ddcx) *2 if for every finitely many bounded 
Borel subsets Bi, . . . , Bn we have that, 

dcx. ddcx) 

It was shown in [4] that verifying the above property for 
all mutually disjoint bounded Borel subsets Bi is a sufficient 
condition for the respective ordering of the pp's. 

Using the above definition, we say that $ is sub-Poisson if 
^ <dcx ^Poi, where ^Poi is some Poisson pp. Noting that for 
the linear function f{x) = x, both / and — / is convex (and 
thus dcx in one dimension) one can observe that $poi needs 
to have the same mean measure as i.e., if $ <dcx ^Poi then 
E[$(B)] = E[$Po,(B)] for every bounded Borel subset B. In 
particular, a stationary pp $ is sub-Poisson if it is dcx smaller 
than Poisson pp ^\ of the same intensity A = E[$(B)]/|B|, 
where is the d —dimensional Lebesgue measure of the set 

B. From now on or will always denote homogeneous 
Poisson pp of intensity A or /i. We will also say that pp 
$ is homogeneous if its mean measure is equal, up to a 



constant, to the Lebesque measure; i.e. E[<1>(_B)] = A|i3|, for 
some constant A all bounded Borel sets B. This is a weaker 
assumption than the stationarity of 

We say that $ is i(ica;(ddca;)-sub-Poisson if $ <idcx{ddcx) 
^Pou where $poi is some Poisson pp. In this case the mean 
measure of $ is smaller or equal (larger or equal) to that of 

^Poi- 



B. Joint Intensities 



The joint intensities '{xi, . . . ,Xk) of a pp are defined by 
the following relation for every finitely many disjoint bounded 
Borel sets Bi, . . . , B„: 



E 



i=l 



■y-Bk 



p^^\xi, . . .,Xk)dxi . ..dxk, 



provided they exist (i.e., the respective moment measures 
admit densities). A pp $ with joint intensities of all orders 
k is said to be weakly sub-Poisson if there exists a constant A 
such that for all fc > 1, 

pW(ari,...,a;fc)<A*^ a.e.. (2.1) 



Due to the idcx property of the function f{xi, . . . ,Xk) — 
Yl - xf, the ordering $1 <idcx ^2 implies that for all > 1, 
E[nti ^i{Bi)] < E[nti «'2(Bi)] for any disjoint bounded 
Borel sets Bi, . . . , Bk. Hence, the respective joint intensities 

(k) 

p) of pp's j — 1,2, provided they exist, obey for 



all k 



> 1, 



p[''\xi, . . . ,Xk) < p'^' {xi, . . . ,Xk) a.e.. Since 
p^"'>{xi, . . . , Xk) = A*^ for Poisson pp $;)^, we are justified in 
using the term weakly sub-Poisson for pp $ satisfying (2.1). 



III. Percolation of Sub-Poisson Gilbert's 
Networks 

Given a pp $ on and a non-negative constant r > 
one defines the (spherical) Boolean model generated by $ of 
ball of radius r as the union C($,r) = Uxe* ^xix), where 
Bx{r) is the ball centered at x e M*^ of radius r. We say that 
C($,r) percolates if there exists an unbounded, connected 
subset of C($, r). This definition extends to any random sub- 
set of W^. One defines the critical radius rc($) of the Boolean 
model as rc($) := inf{r > : P{C($,r) percolates} > 0}. 
Here is the first result of this paper. 

Proposition 3.1: There exist universal constants < c(A) 
and C(A) < 00 depending only on A (and the dimension d) 
such that if <^ <idcx ^\ then c(A) < rc($) and if <ddcx 
thenrd^) < C(A). Thus, for a homogeneous dcx-sub-Poisson 
$ of intensity A we have < c(A) < rc($) < C(A) < 00, 
where the constants depend only on A. The lower bound holds 
also when $ is weakly sub-Poisson. 

The proof is given in the Appendix. The above result can be 
extended to the so-called A: -percolation models. 

Note that the connectivity structure of the Boolean model 
C($,r) corresponds to that of the Gilbert's network with 
nodes in $ and communication range p = 2r; cf. Section 1-Cl. 
Thus the critical communication range pc = Pc(^) in this 
latter network is twice the critical radius of the corresponding 
Boolean model. 

While the finiteness of the critical radius of the Boolean 
model (and thus communication range of the sub-Poisson 



Gilbert's network) is intuitively a desired property, its positiv- 
ity at first glance might be seen from the networking point of 
view as irrelevant, if not a disadvantage. A deeper inspection 
of wireless communication mechanisms shows however that 
sometimes a non-percolation might be also a desired property. 
The following modification of the Gilbert's model, that can be 
seen as a toy version of the SINR model studied in Section IV, 
sheds some light on this latter statement. 

Example 3.2 (Gilbert's carrier-sense network): Consider a 
planar ad-hoc network consisting of nodes modeled by a point 
process $b on the plane (in M^) and having communication 
range p. This process corresponds to a back-bone of the 
network, whose percolation we are looking for. Consider also 
an auxiliary pp $/ of interferers also on M^. Consider the 
following modification of the Gilbert's connection rule: any 
two given nodes of <1>b can directly communicate (are joined 
be an edge) when they are within the communication range 
p from each other, however only when there is no interfering 
node (point of $/) within the sensing range R > p of any 
of these two nodes. Note that any connected component of 
this network is included in some connected subset of the 
complement O = \ C{^i, R — p) of the Boolean model 
C{^j,R — p). Percolation of this vacant region O of the 
spherical Boolean model is thus a necessary condition for 
the percolation of our modification of the Gilbert's network. 
Now, the vacant region of the planar Boolean model cannot 
percolate if the Boolean model itself does (cf [14, Theorem 
4.4]). Thus a non-percolation of some Boolean model related 
to the interferes is a necessary condition for the percolation 
of the communication network on 

IV. Percolation of Sub-Poisson SINR Networks 

In this section, we shall work only on the plane M? . We 
shghtly modify the definition of SINR network introduced 
in [6] allowing for external interferers. The parameters of the 
model are non-negative numbers P (signal power), N (envi- 
rormiental noise), 7 (interference reduction factor), T (SINR 
threshold) and an attenuation function l{r) of the distance 
r > satisfying < l{r) < 1, continuous, strictly decreasing ^ 
on its support, with l{0) > ™ and xl{x)dx < 00. These 
assumptions are exactly as in [7]. 

Given a pp $, the (unit-power) interference generated by 
$ at location x is defined as I^{x) := J2xe^\{x}K\-^ ~ 
x\). More generally this object is also called (additive) shot- 
noise field generated on by $ withe response function I; cf[\, 
Ch. 2.2 and 2.3]. The SINR from x io y with interference 
from $ is defined as 

Pl{\x-V\) 

Let $B and $/ be two pp's on the plane. We do not 
assume any particular dependence between $b and In 
particular one may think of $b C Let P,N,T > 
and 7 > 0. The SINR network with back-bone $b and 
interferers $7 is defined as a graph G($b, ^'7,7) with nodes 
in $B and edges joining any two nodes X, y € <^b when 
SINR(r,X,$/,7) > T and SINR(X, F, $/,7) > T. The 

^So it is rather path gain function. 



SINR(a:,y,$,7) := 



(4.2) 



SNR graph (i.e, the graph without interference) is defined as 
G($b) = G($s,$/,0), which is equivalent also to taking 
= 0. Observe that the SNR graph G'($b) corresponds to 
the Gilbert's network with nodes in $b of communication 
range pi = Percolation in the above graphs is 

existence of an infinite connected component in the graph- 
theoretic sense. 

1 ) Poisson Back-Bone: Firstly, we consider the case when 
the backbone nodes are distributed according to Poisson pp 

^ for some A > 0. We shall use G(A, <!>/, 7) and G(A) 
to denote the corresponding SINR and SNR graphs respec- 
tively. Recall from I-Cl that \c{p) is the critical intensity for 
percolation of the Poisson Gilbert's network of communication 
range p. The following result guarantees the existence of 7 > 
such that for any homogeneous sub-Poisson pp of interferers 
$/ the SINR network G(A, $,7) percolates provided G{X) 
percolates. 

Proposition 4.1: Let A > Xc{pi) and $ <idcx ^nfor some 
p > 0. Then there exists 7c = 'ydX, p, P,T, N) > such that 
G(A, 7) percolates for 7 < 7^. 
The proof is given in the Appendix. 

Recall that we have not assumed the independence of 
and $B = $A- In particular, one can take $/ = ^\U^' where 
$' are some external interferers. If $' is zrfcx-sub-Poisson and 
independent of $a then $/ is also idcx-sub-Poisson (cf [4, 
Proposition 3.2]). The result of Proposition 4.1 in the special 
case of = was proved in [7]. The present extension allows 
for any irfcx-sub-Poisson pattern of independent external in- 
terferers. The proof of our result can be also modified (which 
will not be presented in this version of the paper) to allow for 
external interferer's possibly dependent of the backbone. 
In this fuU generality the result says that any external pattern 
of homogeneous idcx-sub-Poisson interferers added to the 
SINR network of [7] cannot make the giant component of this 
network to disappear, provided the interference cancellation 
factor 7 is appropriately adjusted. Since an irfcx-sub-Poisson 
point process can be of arbitrarily large intensity, the above 
observation can be loosely rephrased in the following form: It 
is not the density of interferers that matters for the network 
connectivity, but their structure; idcx-sub-Poisson interferers 
do not hurt essentially the network connectivity. 

2) Sub-Poisson Back-Bone: We shall now consider the 
case when the backbone nodes are formed by a sub-Poisson 
pp. In this case, we can give a weaker result, namely that 
with appropriately increased signal power P, the SINR graph 
will percolate for small interference parameter 7 > 0. This 
corresponds to an early version of the result for the Poisson 
SINR network, proved in [6], where the percolation of the 
SINR network is guaranteed for the intensity of nodes possibly 
larger than the critical one in the corresponding SNR network. 

Proposition 4.2: Let $b <ddcx for some A > 
and $/ <idcx for some p > and also assume that 
l{x) > for all X > 0. Then there exist P, 7 > such that 
G($B)^/7 7) percolates. 
The proof is given in the Appendix. 

As in Proposition 4.1, we have not assumed the indepen- 
dence of $/ and $B- In particular one can take = $bU$', 
where and $b are independent and dcx-sub-Poisson. 



V. dcx Ordering and Capacity 

Now we want to make some remarks on the relation between 
dcx ordering and some model characteristics usually called ca- 
pacities, including the so-called stochastic-geometric capacity 
functional, as well as some other capacity quantifiers in the 
information-theoretic sense. Our main tool is the following 
result proved in [4, Theorem 2.1]. It says that $i <dcx ^2 
implies the same ordering of the respective shot-noise fields 
(with an arbitrary non-negative response function) 

(/$i(a;i),...,/<i>i(a;„)) <dcx (h^ixi), . . . , I-^^i^n)) (5.3) 

for any Xi,i = 1, . . . ,n. This implies in particular that 



E 



n n 

exp|s^/$i(a;i)|j <E exp^s^U^{xi)^ 



(5.4) 



for both positive and negative s. 



A. dcx Ordering and the Capacity Functional 

Capacity functional Ts{B) of a random set S is defined 
as Th(B) = P{E n B 0} for all bounded Borel 
sets B. A fundamental result of stochastic geometry, called the 
Choquet's theorem (cf [13]) says that the capacity functional 
defines entirely the distribution of a random closed set. The 
complement of it, Vh(B) = 1 — Ts{B), is called void 
probability functional. 

Inequality (5.4) allows to compare capacity functionals and 
void probabilities of Boolean models generated by dcx ordered 
pp's. Indeed, note that 

Vci^,r){B) = e[ [] 1{\X-B\ > r)] = E[exp{ ^ h{X)}] 

where h{x) = log(I(|X - B\ > r)). Note that the latter 
expression has a form of the shot-noise variable and thus 
using the inequality analogous to (5.4) for shot-noise with 
the response function h{-) (cf [4, Theorem 2.1] for such a 
generalization) we observe that if <I>i <dcx ^2 then 



(B) < Vc 
Tc($i,r)(B) > 



C($i,r) 



C(<I.: 



C($2,r) 



){B) 

iB) 



In other words, one can say that dcx smaller pp's exhibit 
smaller voids. In Section VI-B we wiU show some simulation 
examples which illustrate this statement. 

B. dcx Ordering and Network Capacity 

We focus now on capacity quantifiers in the sense of 
communication theory. Inequality (5.4) with s > and 
$1 = $/ and Poisson $2 = will be used in the 
proof of Proposition 4.1 (see Appendix) to show that the 
lower level-sets {x : {x) < M} of the interference field 
generated by $7 percolate through Peierls argument (cf. [1, 
Proposition 14.1.4]) for sufficiently large M. Similarly (5.4) 
with s < can be used to prove that the upper level-sets of the 
interference field 7$^ (x) generated by $/ percolate through 
Peierls argument for sufficiently large level values M. 

Having observed this double impact of sub-Poisson assump- 
tion on $/ it is not evident whether the threshold value 7c of 
the interference reduction factor is larger in sub-Poisson net- 
work than in the corresponding Poisson one. Note that 7c can 
be related to the information-theoretic capacity (throughput) 



that can be sustained on the links of the SINR graph. In what 
follows we will try to explain how sub-Poisson assumption 
can impact some quantifiers of the network capacity. 

1 ) Ordering of Independent Interference Field: It is quite 
natural to consider a network capacity characteristic C = 
/(/(xi), . . . ,/(x„)) that depends on some interference field 
/ through some ddcx function /. We will give a few simple 
examples of such characteristics in what follows. Then, a 
larger in ddcx order interference field (more variable!) I{x) 
leads to larger average capacity E[C]. In particular in the case 
of shot-noise interference field by (5.3) we conclude 
that larger in dcx order pp (clusters more!) leads to larger 
average capacity E[C]. 

Let us illustrate this somewhat surprising observation by 
two simple examples. 

Example 5.1 (Shannon Capacity): Let C(/) = log(l + 
Fq/(N + /)) for some constants Fo,N > and random /. 
Clearly this is a decreasing convex fimction of / and larger 
in convex order / gives larger mean capacity E[C(/)]. The 
smallest value of E[C(/)] given E[/] is attained for constant 

Example 5.2 (Outage capacity of a channel with fading): 
Assume a random channel fading F with convex tail 
distribution function G{t) = P{F > t}. Consider 
C{I) = P{F/{N + I) > T} for some constant 
N,T > 0. Assuming independence of F and / we 
have E[C] = E[G'(r(A^ + /))], which is by our assumption 
expectation of a convex function of / and the same conclusion 
can be made as in Example 5.1. The above general form 
of the expression for the outage capacity can be found in 
many more detailed models based on pp, in particular in 
the Bipolar model of spatial Aloha in [2, Chapter 16]. See 
also [4] for a multidimensional version of this observation. 
Similar conclusion can be made for the ergodic Shannon 
capacity E[log(l + F/{N + /))] of this channel. 

2 ) Ordering of the Back-Bone: Let us take one step further 
and consider the interference created by the original pattern 
of network points (the back-bone process $b) rather than 
an external (independent) interference field. In this case the 
interference (x) at receiver x usually has to be considered 
under the so-called Palm distribution of $5 given the location 
of the emitter (in $b) of this receiver x (again see e.g. the 
Bipolar model of spatial Aloha in [2, Chapter 16] for a detailed 
example). The problem is that dcx ordering of <1>b implies only 
idcx ordering of the respective Palm versions of cf. [4]. 
The fundamental reason for the required "extra" increasing 
property of the comparable functions is that a smaller in dcx 
pp $B (having less clustering or even some point "repulsion") 
will have under Palm probability potentially fewer points in 
vicinity of the conditioned point. This potentially decreases 
interference created locally near this point, thus potentially 
increases our capacity characteristic. More formally: having 
idcx ordered of $5 under Palm probability P° and the 
capacity expressed as a ddcx function of 1^^ we cannot 
conclude any inequality for E°[C(/$g)]. 

The situation is naturally inverted when we are dealing with 
pp which are dcx larger than Poisson pp. A detailed analysis 
in [9] of the outage capacity in the Bipolar model gener- 
ated by some Poisson-Poisson cluster pp known as Neyman- 



Scott pp (which is dcx larger than Poisson pp) confirms the 
above observations. Namely, for smaller transmission distance 
the negative impact of clustering (locally more interferers) 
decreases the outage capacity, while for larger transmission 
distance the positive impact of interferers being more clustered 
increases this capacity. 

3) Multi-hop Capacity Models: Percolation models have 
been also shown useful to study the transport (multi-hop) 
capacity of ad-hoc networks. For example, in [8], by using 
a specific multi-hop transmission strategy that involves per- 
colation theory models, it was shown that a Poisson SINR 
network can achieve capacity rate of the order of l/\/n bits 
per unit of time and per node, thus closing the gap with respect 
to the rate l/^/nlogn shown achievable in Poisson networks 
in [12]. An interesting and open question (particularly in view 
of what was shown above) is whether these capacity results 
can be extended to sub-Poisson networks. 

VI. Examples of Sub-Poisson Point Processes 

From [4, Section 5.2, 5.3], we have a rich source of 
examples of Cox (doubly stochastic Poisson) pp's comparable 
in dcx and ddcx order. In particular, we know that the so 
called Levy-based Cox pp (with Poisson-Poisson cluster pp as 
a special case) is dcx larger than the Poisson pp of the same 
mean intensity. Thus, they can be called as dcx-super-Poisson. 
However, note that any Cox pp, whose (random) realizations 
of the intensity are almost surely bounded by some constant, 
can by coupled with (constructed as a subset of) a Poisson pp 
with intensity equal to this constant. Consequently such Cox 
processes are trivially irfca;-sub-Poisson. 

A. Sub-Poisson Point processes 

In the remaining part of this section we concentrate on the 
construction of examples of dcx-sub-Poisson pp's. 

Example 6.1 (Perturbed lattice): Consider some lattice 
(e.g. the planar hexagonal one, i.e., the usual "honeycomb" 
model often considered in cellular network context). Let us 
"perturb" this "ideal" pattern of points as follows. For each 
point of this lattice, say Zi, let us generate independently, 
from some given distribution, a random number, say iVj, 
of nodes. Moreover, instead of putting these nodes at Xi, 
let us translate each of these nodes independently from Zi, 
by vectors, having some given spatial distribution (say for 
simplicity, of bounded support). The resulting network, called 
perturbed lattice, can be seen as replicating and dispersing 
points from the original lattice. Interesting observations are 
as follows. 

• If the number of replicas Ni are Poisson random vari- 
ables, then the perturbed lattice is a Poisson p.p. 

• If moreover the node displacement is uniform in the 
Voronoi cell of the original lattice, then the resulting 
perturbed lattice is homogeneous Poisson process $a 
with A related to E[A''] and the original lattice density. 

• Now, if Ni are convexly smaller than some Poisson vari- 
able ^ then the perturbed lattice is dcx-sub-Poisson. This 

^i.e., expectations of all convex functions of N = Ni are smaller than the 
respective expectations for Poisson variable of the same mean; a special case 
is when Ni = const is constant 



is SO because the perturbed lattice pp is an independent 
sum of countably many pp formed by the perturbations of 
every vertex of the lattice. This means that ^ = [Jz^^i 
where are independent pp. From [4, Proposition. 3.2], 
we have that if are dcx ordered for every z; , then so 
is Since each of these pp is formed by A^^ i.i.d. pertur- 
bations, it can be shown that these pp's are dcx ordered if 
their respective A^,'s are convexly ordered. This follows 
from proving that the following fiinction g is convex in its 
argument n: g{n) = E(/($2, (^i ),..., (A„))|7V = 
n) for any dcx function /, disjoint bounded Borel subsets 
Ai . . . ,An and any fixed Zi of the original lattice. 
• Consequently, if Ni's are convexly smaller than some 
Poisson random variable and moreover the node displace- 
ment is viniform in the Voronoi cell of the original lattice, 
then the resulting perturbed lattice is rfcx-smaller than the 
respective homogeneous Poisson process ^x. 
The interest in the above perturbed lattice models in the 
networking context stems from the fact that they provide some 
spectrum of models that ranges from periodic grids, usually 
considered in cellular network context, to ad-hoc networks 
almost exclusively considered under Poisson assumptions. In 
Section VI-B we will show some samples of perturbed lattices. 

On the theoretical side, the interest in perturbed lattices 
stems from their relations to zeros of Gaussian analytic func- 
tions (GAFs) (see [15, 16]). More precisely [17] shows that 
zeros of some GAFs have the same distribution as points of 
some "non-independently" perturbed lattice. 

Another class of pp's, which can be shown as weakly 
sub-Poisson, are stationary determinantal pp's. For a quick 
introduction refer [3]. 

Example 6.2 (Determinantal pp's): These pp's are defined 
by their joint intensities satisfying the following relation for 
all > 1 : 

p^^\xi,...,Xk)= Act [k{xi-xM , 

for some Hermitian, non-negative definite, locally square inte- 
grable kernel K ■.M.'^ ^ C Then by Hadamard's inequality 
we have that stationary determinantal pp are weakly sub- 
Poisson. These pp are considered as examples of pp whose 
points "exhibit repulsion". Though zeros of GAFs are related 
to determinantal pp's, curiously enough only zeros of i.i.d. 
Gaussian power series (i.e, f{z) = J2'^=o ^n-^" > with a„ i.i.d. 
standard complex normal) are proved to be determinantal pp's 
(see [15]). 

B. Simulations 

Now, we will show some examples of simulated patterns 
of perturbed lattices. We consider hexagonal pattern of orig- 
inal (unperturbed) points (cf. Figure 1, upper-left plot). The 
replicas are always displaced uniformly in the Voronoi cell 
of the given point of the original point. Consequently, when 
the numbers of replicas N have Poisson distribution then the 
corresponding perturbed lattice is Poisson pp. (cf. Figure 1, 
lower-middle plot). We always take E[A/'] = 1. In order 
to generate sub-Poisson pp's we take N having Binomial 
distribution Bin{n, 1/n). This family of distributions can be 

det{aij)ij < Yli an for Hermitian, non-negative definite matrices (fflij) 



shown to be convexly increasing in n and convex upper- 
bounded by Poisson Poi{l) distribution. Note on Figure 1 that 
in the case Bin{l, 1), where we have only point displacement 
(no replications), the resulting perturbed lattice looks much 
more "regular" (less visible voids, less clusters) than Poisson 
pp. This regularity becomes less evident when n increases and 
already for n > 3 it is difficult to distinguish the perturbed 
lattice patterns from Poisson pp. However, a more detailed 
study of Gilbert graphs on them (not presented here) up to 
approximately n = 5 show that they still have significantly 
different clustering structure than that of Poisson pp. 

The lower-right plot on Figure 1 shows an example of 
perturbed lattice dcx larger than Poisson pp. It is a doubly 
stochastic Poisson pp (Cox pp), where the mean value L of 
the number of repUcas N is first sampled independently for 
each original lattice point, with Pr{ L = a} = 1 — Pr{ L = 
0} = 1/a, for some a > 1 and then, given L, N has 
Poisson distribution Poi{L). This construction is an example 
of the Ising-Poisson cluster pp considered in [4, Section 5.1]. 
Another possibility to generate the number of replicas 
convexly larger than Poisson variable is to "scale-up" the 
Poisson variable of a smaller intensity, i.e., to take A'" = nN' 
where A^' ^ Poi(l/n) for arbitrary n> 1. 

In Figure 2 we consider a larger simulation window com- 
prising about 40^ points of the original hexagonal lattice. 
We generate points of the perturbed lattices with A^ ~ 
Bin{n,l/n) for n = 1,...,3 and n = 20 to approximate 
Poisson pp. For each n we show the Gilbert graph with 
the communication radius p for which the largest component 
(the highlighted one) starts to significantly out-number all 
other components. More precisely for each given simulated 
pattern of points we find the smallest p = p{n) (up to the 
second decimal place) such that the largest component in the 
simulation window contains about 60% of the simulated points 
(cf . the bar-plots showing the empirical fraction of the number 
of points in 10 largest components). Observing values of p{n) 
we conjecture that the critical radius Pc{n) for the percolation 
of the Gilbert graph on the considered perturbed lattices is 
increasing in n. Note that for the unperturbed lattice pc = 1 
(the distance between adjacent nodes in our hexagonal lattice) 
and for Poisson pp of the same intensity pc is known to be 
close to 1.112. A more exhaustive numerical study can be 
found in the extended version of the paper [5]. 

VII. Concluding Remarks 

We have extended two results on nontrivial phase transition 
in percolation of Gilbert's and SINR graphs from Poisson to 
dcx-sub-Poisson pp's. This means that, regarding existence 
of this phase transition, Poisson pp represents a worst-case 
scenario within this class of pp's. 

A natural question in this context is as follows. Consider 
•^i <dcrE ^2- Does this ordering imply that the corresponding 
critical communication ranges in the Gilbert's model are 
ordered as well pd'^i) < Pc(^2)? The answer in the full 
generality might be negative, as we expect some counterexam- 
ples of dcx-super-Poisson pp's with the critical communication 
ranges degenerating to or oo. However, we conjecture that 
the critical communication range of a stationary dcx-sub- 



Poisson pp's of intensity A is smaller than that of a Poisson 
pp of intensity A. 

In this paper we have also discussed the impact of dcx 
ordering on some local quantifiers of the capacity in networks. 
An interesting open question in this context is the relation 
between dcx ordering and the transport capacity in these 
networks. 

Finally, we brought attention to the so called perturbed 
lattices, which can provide a spectrum of pp's monotone in 
dcx order that ranges from periodic grids, usually considered 
in cellular network context, to Poisson pp's and to various 
more clustered pp's including doubly stochastic Poisson pp's. 

Appendix 

In this section we prove our main results of this paper. 

Proof: (of Proposition 3.1) Denote by Z'^(r), the discrete 
graph formed by the vertexes of rZ'^ = {rz : z £ Z'^}, where 
Z'^ is the rf-dimensional integer lattice, and edges between 
z„zj e Z'^ir) such that {z,, + [-^, §]'*}n{z, + [-§, ^]''} ^ 0. 
Now we define a site percolation on Z'^{r) induced by the pp 
$ : Xr{z) = l[^{Qr{z)) > 1] where Qriz) = 
Note that when X2r{z) does not percolate, C($,r) does not 
percolate and if X-!^{z) percolates, so does C($,r). 

The standard technique to show non-percolation is to show 
that the expected number of (self-avoiding) paths of length 
n starting at the origin in the random sub-graph induced by 
opened sites of the percolation model tends to zero as n ^ oo. 
The probability of a path (zi, . . . , z„) of length n in Z^(2r) 
being supported by open sites is Pilli ^(Q2r(-2i)) > 1} < 
E(ni*('32r(-2i))) < (A(2r)'*)" where the first inequaUty is 
due to Markov's inequality and the last inequality is by weak 
sub-Poisson property of $, which in particular is implied by 
/VicT-sub-Poisson property (cf. (2.1)). Since the number of 
paths starting at the origin, of length n in Z'^(2r), is bounded 
by (S"^ — 2)", we have that the expected number of paths of 
length n is at most ((3'' - 2)A(2r)'')" and this tends to 
for r small enough. This shows that there exists c(A) > 
(depending also on the dimension) such that rc($) > c(A). 

For the upper bound, we use the Peierls argument (cf . [11, 
pp. 17-18]) on the site percolation model induced by [z] 
on Z (^). To use this argument, one needs to estimate 
the probability of the site-percolation model not intersecting 
a path (zi,...,z„). This probability can be expressed as 

P{<J>(ULiQ^(^.)) = o}. 

and is smaller than exp{— An(^)''} by the drfca;-sub- 
Poisson property of $ and [4, Proposition 4.1]. Now by 
choosing r large enough this probability can be made as small 
as we wish and so the Peierls argument can be used. ■ 
Proof: (of Proposition 4.1) We follow the proof given 
in [7]. Assuming A > Xc{pi), one observes first that the graph 
G(A) also percolates with any sUghtly larger constant noise 
A'^' = N + 5' , for some 5' > 0. Essential for the original proof 
of the result is to show that the level-set {x : I<j,j (x) < M} of 
the interference field percolates (contains an infinite connected 
component) for sufficiently large M. Suppose that it is true. 
Then taking 7 = 6' /M one has percolation of the level-set 
{y ■ ll^iiy) < S'}. The main difficulty consists in showing 
that G(A) with noise N' = N+5' percolates within an infinite 



hexagonal lattice Bin{l/1) Sm(2, 1/2) Bin{3, 1/3) Poi{l) Cox{5 x Bin{l, 1/5)) 

Fig. 1. "Unperturbed" hexagonal lattice and sub-Poisson perturbed lattices with the number of replicas A'^ having binomial distribution B{n, 1/n). The 
last figure is a super-Poisson perturbed lattice with having double stochastic Poisson (Cox) distribution of random mean L having Bernoulli distribution 
5 X Bm(l, 1/5) (i.e., Pr{ L = 5 } = 1 - Pr{ L = } = 1/5). 
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Fig. 2. Gilbert graph with communication range p and nodes form a perturbed lattice pp with Binomial Bin{n,l/n) number of replicas uniformly 
distributed in hexagonal cells. The largest component in the simulation window is highlighted. Bar-plots show the fraction of nodes in ten largest components. 



connected component of {y : /$j(y) < S'}. This was done 
in [7], by mapping both models G(A) and the level-set of the 
interference to a discrete lattice and showing that both discrete 
approximations not only percolate but actually satisfies a 
stronger, sufficient condition for percolation, related to the 
Peierls argument [1, Proposition 14.1.4]. We follow exactly the 
same steps and the only fact that we have to prove, regarding 
the interference, is that there exists a constant e < 1 such that 
for arbitrary n > 1 and arbitrary choice of locations xi, . . . ,Xn 
one has P{/$j(a;,;) > M, i — l,...,n} < e". To this 
regard, as in [7], using the Chernoff bound we dominate this 
probability by e~''"^^E[exp{^"^j^ s/$^ (xi)}] with arbitrary 
s > 0. The crucial observation for our extension of the 
original proof is that /(ui, . . . , w„) = exp[s(ui + . . . + it„)] 
is an idcx function. By the assumption $/ <idcx and [4, 
Theorem 2.1] we have 



E 



exp|s^/$,(2;i)| <E expjs^ /<j,^^(a;j)| 

i=i " ' 



i=l 



and we can use the explicit form of the Laplace transform 
of the Poisson shot-noise (the right-hand-side in the above 
inequality) to prove, exactly as in [7], that for sufficiently small 
s it is not larger than K"" for some constant K which depends 
on p but not on M. This completes the proof. ■ 

Proof: {of Proposition 4.2) In this scenario, increased 
power is equivalent to increased radius in the Gilbert's model 
associated with the SINR model. From this observation, it 
follows that by using the discrete mapping and arguments as 
in the proof of Proposition 3.1, we obtain that with increased 
power the associated Gilbert's model percolates. Then, we use 
the approach from the proof of Proposition 4.1 to obtain a 
7 > such that the SINR network percolates as well. ■ 



References 

[1] F. Baccelli and B. Blaszczyszyn. Stochastic Geometry and Wireless 

Networks, Volume I — Theory, volume 3, No 3^ of Foundations and 

Trends in Networking. NoW Publishers, 2009. 
[2] F. Baccelli and B. Blaszczyszyn. Stochastic Geometry and Wireless 

Networks, Volume II — Applications, volume 4, No 1-2 of Foundations 

and Trends in Networking. NoW Publishers, 2009. 
[3] J. Ben Hough, M. Krishnapur, Y. Peres, and B. Virag. Determinantal 

processes and independence. Probability Surveys, 3:206-229, 2006. 
[4] B. Blaszczyszyn and D. Yogeshwaran. Directionally convex ordering of 

random measures, shot-noise fields and some applications to wireless 

networks. Adv. Appl. Probab., 41:623-646, 2009. 
[5] B. Blaszczyszyn and D. Yogeshwai'an. Connectivity in sub-Poisson 

networks, http://hal.inria.fr/inria-00497707, 2010. 
[6] O Dousse. F. Baccelli, and P Thiran. Impact of interferences on 

connectivity in ad-hoc networks. IEEE/ACM Trans. Networking, 13: 

425-543, 2005. 

[7] O. Dousse, M. Franceschetti, N. Macris, R. Meester, and P. Thiran. 
Percolation in the signal to interference ratio graph. Journal of Applied 
Probability, 43(2):552-562, 2006. 
[8] M. Franceschetti, O. Dousse, Tse D. N. C, and Thiran P. Closing the 
gap in the capacity of wireless networks via percolation theory. IEEE 
Trans. Inf Theory, 53:1009-1018, 2007. 
[9] R. Ganti and M. Haenggi. Interference and outage in clustered wireless 
ad hoc networks. IEEE Tr Inf Theory, 55:4067^086, 2009. 
[10] E. N. Gilbert. Random plane networks. SI AM J., 9:533-543, 1961. 
[11] G. R. Grimmett. Percolation. Springer- Verlag, New York, 1989. 
[12] P. Gupta and P. R. Kumar The capacity of wireless networks. IEEE 

Transactions on Information Theory, 46(2):388^04, 2000. 
[13] G. Matheron. Random Sets and Integral Geometry. John Wiley & Sons, 
London, 1975. 

[14] R. Meester and R. Roy. Continuum percolation. Cambridge University 
Press, 1996. 

[15] Y. Peres and B. Virag. Zeros of the i.i.d. Gaussian power series: a 
conformally invariant determinantal process. Acta Mathematica, 194: 
1-35, 2005. 

[16] M. Sodin and B. Tsirelson. Random complex zeroes; I. asymptotic 

normality. Israel J. Math., 144:125-149, 2004. 
[17] M. Sodin and B. Tsirelson. Random complex zeroes; II. perturbed 

lattice. IsraelJ. Math., 152:105-124, 2006. 



